Very high precision bound state spectroscopy near a 85 Rb Feshbach resonance 
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We precisely measured the binding energy (tbind) of a molecular state near the Feshbach resonance 
in a 85 Rb Bose-Einstein condensate (BEC). Rapid magnetic field pulses induced coherent atom- 
molecule oscillations in the BEC. We measured the oscillation frequency as a function of B-field and 
ht the data to a coupled-channels model. Our analysis constrained the Feshbach resonance position 
[155.041(18) G], width [10.71(2) G], and background scattering length [-443(3) a ] and yielded new 
values for vds, vdt, and C§. These results improved our estimate for the stability condition of an 
attractive BEC. We also found evidence for a mean-held shift to ebind- 
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The phenomenon of a Feshbach resonance in ultracold 
collisions of alkali atoms has received much theoretical 
and experimental interest in recent years and has sparked 
interest in the subject of resonant Bose-Einstein conden- 
sates (BEC). Feshbach resonances 0,0 have been used to 
control elastic and inelastic collisions in ultracold gases 

g la, lol and for tuning the self-interaction in BEC 
by changing the magnitude of an exter- 
nal magnetic field. The magnetic field controls the self- 
interaction in the BEC by affecting the s-wave scattering 
length, a. Close to resonance, the scattering length varies 
with B-ficld according to 
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where B poa k is the resonance position and is defined to 
be the magnetic field where the magnitude of a becomes 
infinite, a Dg is the background scattering length, and 
A = B zeTO — Speak is the resonance width where B zcro is 
the B-field where the scattering length crosses zero. Mea- 
surements of Feshbach resonance positions and widths 
have been used in a variety of alkali atoms to improve the 
determination of the interatomic potentials. These po- 
tentials have then been used to precisely calculate a mul- 
titude of im port ant properties for trapped atomic gases 

Recently, we applied rapid magnetic field variations 
near a Feshbach resonance to create an atom-molecule 
superposition state in a 85 Rb BEC ^(|, which has al- 
lowed us to precisely determine the Feshbach resonance 
position and width. Our novel technique for studying 
the Feshbach resonance relies on the presence of atom- 
molecule coherence 0, 0, 0. By inducing periodic 
oscillations in the number of condensate atoms, we ob- 
tain a direct, high precision measurement of the molec- 
ular bound state energy. Exploiting the resonance, we 
tune the molecular state very close to threshold — to 
our knowledge, this is the most weakly bound state ever 
observed. The present method for studying the Feshbach 
resonance through atom-molecule oscillations offers all of 



the many inherent advantages of a frequency measure- 
ment, including the possibility of high measurement pre- 
cision, a lack of sensitivity to errors in the absolute atom 
number calibration, and a simple interpretation of the 
oscillation frequency in terms of the relative energy dif- 
ference between the atomic and molecular states. When 
these advantages are combined with an improved method 
for magnetic field calibration |2(|, the present technique 
for probing the Feshbach resonance is much more precise 
than previous experiments that examined such Feshbach 
resonance obscrvables as variable rethcrmalization rates 
in a trapped cloud of atoms 0,0], enhancements of pho- 
toassociation rates and inelastic loss rates near the 
resonance |2lj . and variations of the mean-field expan- 
sion energy of a BEC • 

To complete our precise characterization of the Fesh- 
bach resonance, we also made an improved measurement 
of B zero , the magnetic field where the scattering length 
vanishes. This experiment is very similar to our previ- 
ous work [Til |2^ . where we determined the a=0 field 
by measuring the critical number (N cr ;t) for collapse of 
a BEC, and then we extrapolated to the magnetic field 
where N cr ; t would be infinite. We have improved the mea- 
surement precision by about a factor of 4 by improving 
our magnetic field calibration and using a larger number 
of condensate atoms to measure the collapse. We find 
B zoro =165. 750(13) G. 

The procedure used to generate atom-molecule os- 
cillations in 85 Rb Bose-Einstein condensates has been 
described in previous work |l6j |. so we merely outline 
the method here. After creating condensates with ini- 
tial number of atoms No ~ 16000 at a magnetic field 
B~162 G, we apply two short B-field pulses (~40 fis du- 
ration) that approach and then recede from the Feshbach 
resonance at Bp^ ~155 G. The intermediate value of 
magnetic field between the pulses, B ovo i vc , and the time 
spacing between pulses, t OV oivc, are variable quantities. 
The double pulse sequence is followed by a slow change 
in the B-field to expand the BEC @, then the trap is 
switched off (B— >0) and destructive absorption imaging 
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is used to count the number of atoms remaining in the 
condensate. 

As in Ref . 01 > periodic oscillations in the BEC number 
were observed as a function of t ovolvo (see Fig. 1). We 
fit the BEC number oscillation to a damped harmonic 
oscillator function with an additional linear loss term: 

N(t) = N avg - at + Acxp (-0t) sin (u> e t + c/>), (2) 

where N avg is the average number, A is the oscillation 
amplitude, a and /3 are the number loss and damping 
rates, respectively, and u> e = 2-k^Jvq — [j3/(2ir)] 2 . The 
quantity of interest here is v$, the natural oscillator fre- 
quency corresponding to the molecular binding energy, 
"o =e bind/ h- We measured the oscillation frequency for 
values of B cvo i vc from 156.1 G to 161.8 G. Over this 
range, the frequency varies by over 2 orders of magni- 
tude (10-1000 kHz), but the linear loss rate changes very 
little. The damping rate shows a significant B-field de- 
pendence, increasing from /3 ~ 2ir x 0.8 kHz near 156 G to 
(3 ~ 2tt x 22 kHz near 162 G. We find no significant den- 
sity dependence to the damping at B ovo i vo =158.60(5) G; 
increasing the total atom density by a factor of 4.3(5) 
leads to a damping rate increase of only a factor of 1.3(3). 
Atom loss from the BEC is well described by a linear rate 
of -2 to -7 atoms//xs over the field range of interest. The 
rate is consistent with previous measurements of number 
loss due to a single B-ficld pulse toward the Feshbach 
resonance |2c| . 

To characterize the Feshbach resonance, it is neces- 
sary to know both the oscillation frequency and B evo i ve . 
We precisely measured B evo i ve by transferring atoms to 
an untrapped spin state by driving Am=+1 spin flip 
transitions with an applied pulse of rf radiation (pulse 
lcngth=5^25 /is). The spin flip frequency was deter- 
mined from the rf lincshapc for the loss of atoms from 
the magnetic trap. After measuring the rf transition fre- 
quency, we inverted the Breit-Rabi equation to obtain 
the corresponding B-ficld. To ensure that the magnetic 
field was sufficiently constant during t cvo i ve , we mapped 
out B(t) using rf pulses with lengths short compared to 
tovolvo- Due to interference of the rf radiation with the 
magnetic field control circuitry, there was a small sys- 
tematic shift of the field as a function of the rf power 
used. The uncertainty for each magnetic field determi- 
nation was the quadrature sum of the uncertainty due 
to the lincshapc measurements (~15 inG) and the uncer- 
tainty in the extrapolation to zero rf power (—20 mG). 
Typically, the total uncertainty in the average B-field was 
-25 mG. 

The measured oscillation frequencies versus magnetic 
field are plotted in Fig. 2. We use these data and the 
zero-crossing field B zcro to completely characterize the 
scattering length and binding energy as a function of 
magnetic field near the Feshbach resonance. As a starting 
point we use the coupled-channels analysis of van Kem- 
pen et al. 0], where several high-precision data for 85 Rb 
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FIG. 1: BEC number versus pulse spacing, t cvo i vc . (a) 
B cvo i vc =156. 840(25) G. At this magnetic field, which is 
relatively close to resonance, the oscillation frequency is 
very low (i/o=9. 77(12) kHz) so that the damping and 
atom loss significantly affect the observed time dependence 
(here /3=2tt x 0.58(12) kHz and a=7.9(4) atoms/^s). (b) 
B ovo ivc=159. 527(19) G. Farther from resonance, the time de- 
pendence of the BEC number is dominated by the higher 
oscillation frequency of ^o = 157.8(17) kHz. Damping of the 
oscillations and atom loss are negligible in the relatively short 
time window used to determine vq. 



and 87 Rb were combined to perform an inter-isotope de- 
termination of the rubidium interactions with unprece- 
dented accuracy. The predictive power of this analysis 
can be seen from Ref. |l6j , where the initial data on the 
atom-molecule coherence were already in good agreement 
with the predicted binding energy of the underlying Fes- 
hbach state. Another example of the accuracy of the 
analysis in Ref. is its agreement with more than 40 
Feshbach resonances recently discovered in 87 Rb pjj . 

Van Kempen et al. used the best known values Q 
for the resonant magnetic field -B poa k and zero crossing 
B ZC ro- In this Letter, we ignore the relatively imprecise 
value of Bpoak from Ref. |hj| , and instead use the mea- 
sured dependence of binding energy on magnetic field 
along with the new B zero measurement given above to 
determine the interaction parameters. We observe that 
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FIG. 2: Molecular binding energy versus magnetic field, 
Bcvoivo. (a) The points are measured values of the atom- 
molecule oscillation frequency, vo, while the solid line repre- 
sents the molecular binding energy, which we fit to the data 
by adjusting the parameters of a coupled-channels scatter- 
ing theory. Only black points were included in the fit; white 
points were excluded because they experienced a statistically 
significant mean-field shift. To improve visibility, the points 
are larger than the error bars, (b) Same as in (a), but with 
a linear scale for the vertical axis. The inset shows the devi- 
ation of the lowest frequency data from the fit to the rest of 
the data. 



the fitting procedure is mainly sensitive to only three pa- 
rameters: the van der Waals dispersion coefficient, Cq, 
and the non-integral vibrational quantum numbers at 
dissociation, vds and vdt, which determine the position 
of the last bound state in the singlet and triplet poten- 
tials, respectively. Varying the additional parameters Cs, 
Cio , (ft j, (the first-order energy-dependence of the phase 
of the oscillating triplet radial wave function), and J, 
the strength of the exchange interaction, does not im- 
prove the fitting because these changes can be absorbed 
in small shifts of vus, vdt, and Cq. Therefore, we take 
the mean values for these four parameters [24| from the 
most recent determination in Ref. |23j . 

The best fit to B zero and the seven highest frequency 
data points yields a reduced x 2 =0.30 for 5 degrees of 



TABLE I: Sensitivities of the determined interaction param- 
eters vds, vdt and Cs to fractional uncertainties in Cg, Cio, 
<p® and J. For instance, the systematic error in C% due to a 
10% uncertainty in Cg is 123 x 0.10 = 12.3 a.u. 
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freedom. This value of \ 2 is improbably low due to 
the fact that the uncertainty in the data is dominated 
by the systematic uncertainty in magnetic field related 
to the magnitude of the rf power shift. Figure 2 shows 
the theoretical fit to the binding energy data as a func- 
tion of magnetic field. From the fit, we find substan- 
tially improved values for the Feshbach resonance po- 
sition B peak =155.041(18) G, width A=10.71(2) G, and 
background scattering length ab g =-443(3) a . These 
results may be compared to previously obtained re- 
sults B pcak =154.9tt) G and A=11.0(4) G Q, and 
a bg=-450(3) ao |25j. Our best interaction parame- 
ter values are C 6 =4707(2) a.u., Vr,g=0.00918(17), and 
vdt=Q- 94659(29). Here the error bars do not include 
systematic errors due to the uncertainties in other inter- 
action parameters that are not constraine d by our data. 
To compare our values with those of Ref. |l5j , we deter- 
mined the sensitivity of our three interaction parameters 
to systematic shifts in the other parameters, as shown in 
Table I. Using the fractional uncertainties in Cs, Cio, 
</>|, and J from Ref. 0, we find C 6 =4707(13) a.u., 
ij_ds=0.0092(4), and u DT =0.9466(5V All of these val- 
ues agree with those given in Ref. |l5j: C6=4703(9) a.u., 
V£>s=0.009(l), and v flT =0.9471(2_L Our value for v DS 
is more precise than that of Ref. [l5j, while vdt and Cq 
arc slightly less precise. If future experiments allow im- 
provements to the other interaction parameters, then our 
results will also become more precise since the system- 
atic errors are comparable to or larger than our statistical 
errors from the fit. 

To understand the strong parameter constraints that 
we obtain with our bound state spectroscopy, it is impor- 
tant to consider the nonlinear dependence of the binding 
energy on magnetic field. The magnetic field dependence 
of ebind as it approaches the collision threshold depends 
sensitively on the exact shape of the long range inter- 
atomic potentials, which arc mainly characterized by the 
van der Waals coefficient, Cq. At magnetic fields far from 
resonance, the bound state wave function is confined to 
short internuclear distance and the binding energy varies 
linearly with magnetic field. The linear dependence on 
B-field gives relatively little information about Cq. As 
the B-field approaches resonance, the detuning decreases 
until the bound state lies just below threshold. Now the 
bound state wave function penetrates much deeper into 
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the classically forbidden region, which causes ebmd to 
curve toward threshold as a function of magnetic field. 
Because the energetically forbidden region stretches out 
as Ce/r 6 , the observed curvature depends sensitively on 
the Cq coefficient. One can show |2(j that an analytical 
Feshbach model that includes the correct potential range 
and background scattering processes can reproduce 
the binding energy curve over the full range of magnetic 
field. 

The coupled-channels theory used in this work applies 
to two-body scattering; therefore, this theory cannot ac- 
count for many-body effects in the atom-molecule BEC 
system, such as a mean-field shift to the observed oscilla- 
tion frequency 0,H3|- Any such mean-field shift must be 
fractionally largest near the Feshbach resonance, where 
the binding energy approaches zero while the atom-atom 
scattering length increases to infinity. We searched for a 
mean-field shift to the oscillation frequency when B evo i ve 
was decreased to ^156 G. As shown in Fig. 2, the low- 
est magnetic field data display a clear frequency shift 
with respect to the coupled-channels theory prediction. 
As Bovoivo approaches resonance, the observed shift in- 
creases to 1.7 kHz, which significantly exceeds a simple 
estimate for the average atom-atom mean-field shift in 
the BEC: ^h 2 (n)a/m ~ 0.5 kHz at B cvolvc =156.1 G. 
We are presently investigating new experimental tech- 
niques to further study the mean-field shifts, including 
their density dependence. 

Since the lowest frequency data show evidence for a 
mean- field shift, we exclude these points from the (two- 
body) theory fit. We determine the cutoff magnetic field 
for the excluded region by the following procedure. We 
fit the data set that includes all frequency measurements 
satisfying isq > 9 kHz. Eliminating the lowest frequency 
point from the set causes the reduced x 2 to decrease from 
0.3 to 0.2, and there is no significant change in parame- 
ter values. In contrast, adding the next lower frequency 
point increases the reduced x 2 to 1.9, causing a system- 
atic shift in the parameter values. The observed behav- 
ior seems sensible since we expect mean-field shifts to in- 
crease rapidly as one moves toward resonance (see Fig. 2). 

As a result of the improved determination of the 
85 Rb Feshbach resonance parameters, we find that our 
new value for the off-resonant or background scattering 
length, <2b g =-443(3) ao, is inconsistent with the value 
given in Rcf. |14|. where ab g =-380(21) a . The most 
plausible explanation we can find for disagreement is that 
the theoretical expression used to relate measured rether- 
malization rate to cross section is insufficient for the req- 
uisite level of accuracy. However, the new value for ab g 
allows us to revise our previous estimate for the stability 
condition of a BEC with negative scattering length [22|. 
We use Eq. to obtain the linear slope of scattering 
length versus B-field near B=B zoro . We then find the 
stability coefficient for BEC collapse, fc C oiiapse, by com- 
bining the value of Aa/AB=-39. 87(22) a /G with the 



measured slope of 1/N cr ;t versus magnetic field [22J of 
0.00126(3) (atoms G) -1 . Thus, we obtain the revised 
value fccoiiapsc^O. 547(58), where the error is dominated 
by a 10% systematic uncertainty in the determination of 
N cr j t . The present determination agrees with the theo- 
retical value of 0.55 [29j . 

In conclusion, we present a unique method for explor- 
ing a 85 Rb Feshbach resonance. The observed atom- 
molecule coherence allows us to study the highly nonlin- 
ear dependence of the molecular binding energy on mag- 
netic field. We find good agreement with an analysis of 
van Kempen et al. [15j and improved precision for the 
characterization of the Feshbach resonance. In addition, 
we observe mean-field shifts to the binding energy, offer- 
ing the possibility for future studies of many-body effects 
in this exciting system. 
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